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Abstract 

We consider the following problem: When do alternate eigenvalues taken from a matrix ensemble them- 
selves form a matrix ensemble? More precisely, we classify all weight functions for which alternate eigenvalues 
from the corresponding orthogonal ensemble form a symplectic ensemble, and similarly classify those weights 
for which alternate eigenvalues from a union of two orthogonal ensembles forms a unitary ensemble. Also 
considered are the fc-point distributions for the decimated orthogonal ensembles. 

1 Introduction 

Given a probability measure on a space of matrices, the eigenvalue PDF (probability density function) follows 
by a change of variables. For example, consider the space of n x n real symmetric matrices A = [flj fc]o<j fc<„ 
with probability measure proportional to 



e "-{dA), (dA) :=Y[a,k (1.1) 



-Tr(A^)/2 _ 

The eigenvalues xq < ■ ■ ■ < Xn-i are introduced via the spectral decomposition A — RLEF where R is 
a real orthogonal matrix with columns given by the eigenvectors of A and L = diag(xo, . . . , x„_i). Since 
g-Tr(A )/2 _ g-Ej=o ^jjg change of variables is immediate for the weight function; however the change of 
variables in {dA) cannot be carried out with such expedience. 

The essential point of the latter task is to compute the Jacobian for the change of variables from the 
independent elements of A to the eigenvalues and the independent variables associated with the eigenvectors. 
Also, because only the eigenvalue PDF is being computed, one must integrate out the eigenvector dependence. 
In fact the dependence in the Jacobian on the eigenvalues separates from the dependence on the eigenvectors. 
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so the task of performing the integration does not become an issue. ExpUcitly, one finds (see e.g. pQ|) 



(dA) = A{x) Jl dxj{R^dR) 



where A(a;) := A(xo, . 
proportional to 



0<j<k< 



), and thus the eigenvalue PDF corresponding to (|]^) is 



l[g{x,)\A{x)\ 



(1.2) 



(taking the absolute value of A(x) allows the ordering restriction on the eigenvalues to be 



with g{x) 
dropped). 

More generally, the above working shows that a space of n x n real symmetric matrices with probability 
measure proportional to 



exp I 



(1.3) 



will have eigenvalue PDF (F2) with g{x) = exp(^°^-^ ajX^). Because (O) is unchanged by similarity trans- 
formations A h-> RAR^ with R real orthogonal, for general g (1.2) is said to be the eigenvalue PDF of an 
orthogonal ensemble, and denoted OE„(g). 

Real symmetric matrices are Hermitian matrices with all elements constrained to be real. If one considers 
Hermitian matrices without this constraint, so the off diagonal elements can now be complex, the eigenvalue 
PDF corresponding to the ensemble (1.3) is again given by (1.2) but with |A(a;)| replaced by (A(a;))^. Because 
(L3) is then unchanged by A i-^ UAU'^ for U unitary, (L2) so modified is referred to as a unitary ensemble 
and denoted UE„(5). The third and final possibility |^] is to consider n x n Hermitian matrices in which each 
element is itself a 2 x 2 matrix of the form 



z w 
-w z 



(1.4) 



This class of 2 x 2 matrices form the real quaternion number field H. The spectrum of such matrices, regarded 



as 2n X 2n matrices with complex entries, is doubly degenerate. The ensemble of matrices (1.3) is now invariant 
under the transformations A i-^ BAB'' for B symplectic unitary, and so referred to as a symplectic ensemble. 
The eigenvalue PDF of the distinct eigenvalues is given by (1.2) with |A(a:)| replaced by (A(a;))^, and this is 
denoted SE„(5). 

The matrix ensembles corresponding to the eigenvalue PDFs 

0E„(e--'/2), UE„(e-^'), SE„(e--') 

are given the special labels GOE„, GUE„ and GSE„ respectively (the G standing for Gaussian). As seen from 
(1.1) they can be realized by an appropriate Gaussian weight function in the probability space. Because for A 
real symmetric 



n— 1 n— 1 

-Tr(A^)/2 ^ 77 ,-4, -Q , 

j<k 
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and similarly for A Hermitian with complex or real quaternion elements, independent elements of the Gaussian 
ensemble are independently distributed Gaussian random variables. 

There are also a number of other known random matrix ensembles with this latter property, and which have 
eigenvalue PDF of the form 0E„(5), \JEn{g) or SE„(5) for some g. Seven such ensembles result by taking the 
Hermitian part of the matrix Lie algebras related to Cartan's ten families of infinite symmetric spaces [|9) . We 
specify five of these: 



Mat(p,g;: 
Mat(p, q; I 
Mat(p, q;] 
Symm(n; ' 
Anti(n; ' 



p X q matrices over M {p > q) 

p y. q matrices over C {p > q) 

p X q matrices over H {p > q) 

n X n symmetric complex matrices 

n X n antisymmetric complex matrices. 



(1.5) 
(1.6) 
(1.7) 
(1.8) 
(1.9) 



The quantities of interest are the square of the non-zero singular values, or equivalently the eigenvalues of 
A for A a member of the ensemble, in each case. The first two of these ensembles were studied long ago in 
mathematical statistics [^8[ |l^ ; these two together with the third have occured in recent physical applications 
(see Q and references therein), while the final two (in a different guise) have also arisen in a physical context 
The distribution of the eigenvalues of A'^ A can be computed in a number of ways; one approach is to make 
use of the correspondence to a symmetric space (of types BDI, AIII, CII, CI and Dili respectively), 
which allows the tables in [|l5| to be utilized. Abusing notation, we have 



Mat(p,(7;R 
Mat(p, q; C 
Mat(p, q; H 
Symm(n; C 
Anti(2n; C 
Anti(2n + 1;C 



OE,(a; 



2{p-q) + l-x 



= 0E„(e--/2) 
= SE„(e--) 



(1.10) 



Up to the scale of x, all the above weight functions are of the Laguerre form a;"e~^ and so by definition are 
examples of Laguerre matrix ensembles. 

Another class of matrix ensembles in which the entries of the underlying matrices are independently dis- 
tributed Gaussian random variables are known in mathematical statistics |2^. With a <E Mat(pi, q; F) (where 
F = R, C or H), b E Mat(p2, q,^), and A = a^a, B = b^b, these distributions are described by 

Beta(pi,p2, <z; IF) qx q matrices A{A + B)^^ . 

They have corresponding eigenvalue PDF (abusing notation as in ( |l.lO| )) 

Beta(pi,p2,g;K) = 0E,(a;(Pi-«-i)/2(^ _ ^)(P2-g-i)/2) 

Beta(pi,p2,g;C) = UE,(xPi"«(l 



Beta(pi,p2,'?;IHI) = SEg{x' 



2{pi-q) + l 



(l-x) 



2ip2-q) + l 



(1.11) 
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where < a: < 1, and thus involve weight functions of the Jacobi type. 

The above revision demonstrates that it is possible to realize, in terms of matrices with entries which are 
independently distributed Gaussian random variables, the distributions OE.,i(g), \JF,n{g) and SF,n{g) for g one 
of the forms 

e""', a;"e-", x^il-xf. (1.12) 

These same weight functions occur in the theory of orthogonal polynomials |^ — they are associated with the 
three families of classical orthogonal polynomials Hermite, Laguerre and Jacobi respectively, and are themselves 
referred to as classical weight functions. The classical polynomials share many special properties not enjoyed 
by orthogonal polynomials associated with other weight functions. In the present study of matrix ensembles, 
we will see that the distributions OE„(g), UE„(g) and SE„(f/) also have special features for g a classical weight 
function (1.12). 

Our interest is in the properties of alternate eigenvalues in matrix ensembles. In particular we seek to 
determine the weights g for which alternate eigenvalues taken from a random union of two orthogonal ensembles 
form a unitary ensemble. Similarly we seek the weights g for which alternate eigenvalues from an orthogonal 
ensemble form a symplectic ensemble. The motivation for this study comes from recent work of Baik and Rains 
[Q. Consider the distribution OE„(e~'^), n even, and order the eigenvalues xq < xi < ■ ■ ■ < a;„_i. In Q it 
was proved that after integrating out every second eigenvalue Xn-i,Xn-3, ■ ■ ■ etc. the remaining eigenvalues 
have the distribution SE„/2- The proof of Baik and Rains is particular to the a = case of the Laguerre 
ensemble. However other considerations lead these authors ^ to conjecture that in an appropriate scaled limit 
the distribution of the largest eigenvalue in the GSE corresponds to that of the second largest eigenvalue in 
the GOE. From this it is remarked that presumably the joint distribution of every second eigenvalue in the 
GOE coincides with the joint distribution of all the eigenvalues in the GSE, with an appropriate number of 
eigenvalues. 

Baik and Rains |^ were also led to consider two GOE„ spectra, superimposing them at random, and 
integrating out every second eigenvalue of the resulting sequence. Results were presented which suggest that in 
the scaled n ^ oo limit at the soft edge the distribution becomes that of GUEoo, appropriately scaled. 

Such inter-relationships between ensembles first occured in the work of Dyson on the circular ensembles 
of random unitary matrices. This ensemble has eigenvalue PDF proportional to 

Y[ |e'«'= -e^^^l'', 0<9^<2n (1.13) 

0<j<fe<n 

for P — 1,2 and 4 (COE„, CUE„ and CSE„) respectively. Dyson conjectured that 

alt(COE„ U COE„) = CUE„ (1.14) 

which means that if two spectra from the COE„ distribution are superimposed at random with every second 
eigenvalue integrated out, the CUE„ distribution results. This was subsequently proved by Gunson |l^. Also, 
Mehta and Dyson |19 proved that integrating out every second eigenvalue from the distribution COE„ with n 
even gives the distribution CSE„/2i or symbolically 

ah(COE„) = CSE„/2 . (1.15) 
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The circular ensembles can be analyzed in the course of the present study of ensembles with real valued 
eigenvalues by making the stereographic projection 



1 + iXA 



The PDF (1.13) then maps to 



Ha 



1 



(l + a;2)'3("-i)/2+i 



n W-^3f 



j=0 ^ ■ ' J ' Q<j<k<n 

which is of the general type under consideration. Here the weight function is of the form 

1 

a> 1. 



(l+x2)"^ 



(1.16) 



This only has a finite number of well defined moments and thus in this respect differs from the classical weight 
functions (1.12). On the other hand the corresponding orthogonal polynomials are °'\ix)}^ 



with pIC'^^ denoting the Jacobi polynomial, thus implying (1.16) can be viewed as a fourth classical weight 
function. 



l<Q-l/2 



2 Pseudo-ensembles 

We begin with the orthogonal ensemble eigenvalue PDF ( |l.2| ), taking away the modulus sign, replacing n by I 
(to avoid overuse of the former) and rewriting the product as a determinant using the Vandermonde formula 
to obtain 

^(2;) ]^.9(a^i) = det(x])o<jj<i ]^g(xi) = det(g(xi)a;^)o<ij<i. (2.1) 

i i 

In particular, we note that each row corresponds to a variable, while each column corresponds to a function. 
Given a collection of n functions i^; : R ^ K, we thus define the associated "orthogonal pseudo-ensemble" by 
the following "density" : 

det(^;(a;,))o<»j<i- (2.2) 

Thus any orthogonal ensemble is also an orthogonal pseudo-ensemble, but certainly not vice versa. Indeed, one 
has: 

Theorem 2.1. Fix an integer I > Q, and Zet G : M ^ R 6e a function supported on at least n points. Then for 
a collection of n functions Fq, Fi, . . . -Fi-i, we have 

det(Fj(xj))o<ij<i oc []G(xOA(x) (2.3) 

i 

if and only if there exist I linearly independent polynomials pi of degree at most l~ 1 such that Fi{x) — pi{x)G{x) 
for all i and x. 
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Proof. The "if" portion is easy enough: 

det{G{xi)pj{xi))a<ij<i = ]^ G(a;i) det(pj(xi))o<j,j</ « G(xi)A(a;), (2.4) 

i i 

since the polynomials are assumed linearly independent. 

Now, suppose (2.3) holds. It will turn out to be convenient to restate the equation in terms of exterior 
products. Define a vector- valued function Vf{x) by 

Vf{x),=F,{x). (2.5) 

Then we can write 

det(F,(x,))o<^,,<i - (V>(xo),( /\ Vpixi))), (2.6) 

l<i<i 

where (, ) stands for the standard duality between 1-forms and I — 1-forms. Consider this as a function of xo as 
the other variables range over the support of G; we have: 

Vf{x) ■ ( f\ Vf{x,)) cx G{x) \{{x- X,). (2.7) 
i<i<i i<i<i 

Now, since G has at least I elements in its support, these functions span an Z-dimensional space (this follows, 
for instance, from Lagrange's interpolation formula). On the other hand, the functions must clearly be linear 
combinations of the Fi . Since there only I functions Fi , it follows that we can write the Fi as linear combinations 
of the functions G{x)x^, Q < i < I. But this is precisely what we wanted to prove. □ 

Similarly, the density function of a symplectic ensemble can also be written as a determinant, namely 

^(2;)''n5(^*)^ = det(5(a;j)a;^,j5(xj)a;^~^)o<j<Lo<j<2i; (2.8) 

i 

this follows by differentiating the Vandermonde determinant. When log{g) is differentiable, we can perform 
column transformations to put this determinant in the form 

det{Fj{x,),F-{x,))o<t<LQ<j<2i; (2.9) 
simply take Fj{x) = g{x)x^ , and observe that 

F;ix)^^F,{x)^jgix)x^-K (2.10) 

In fact, we can often define even when the functions Fj are not differentiable, by expressing it in terms of 
the 2-form-valued function 

Vj,^\x) = lim -^—(VFix) A VFiy)). (2.11) 

y^x l^x ~ y) 

For Fj(x) = g(x)x^ , we find that this is defined wherever g is continuous. 
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Theorem 2.2. Fix an integer I > 1, let O be a nonempty open subset o/M, and let G : O M be a continuous 
function supported on O. Then for a collection of continuous functions Fj : O ^ M. such that ( ^.9[ ) is well- 
defined, 

A(x)4J]G(xO' (2.12) 

i 

on O^' if and only if there exist linearly independent polynomials pj of degree at most 21 — 1 with Fj{x) = 
G{x)pj{x). 

Proof. Again the "if" case is straightforward. In the other direction, we can clearly divide each Fj by G, and 
thus may assume WLOG that G = 1 on O. 
We first consider the case 1 = 2, for which 

dct{F,{x) F;{x) F,{y) F;{y))o<,<i = {vj.'\x), vj.^\y)) = [x - yf. (2.13) 

As y varies over O, this spans a 5-dimensional function space; it follows that as y varies, Vp (y) spans a 
5-dimensional space (the dimension must be either 5 or 6; 6 clearly leads to a contradiction). In other words, 
there must be a linear dependence between the coefficients of Vp {y). By replacing the Fi with an orthogonal 
linear combination, we find that this dependence is WLOG of the form 

^i'^(2/)oi=Gyi')(y)23, (2.14) 

for some constant G. Now, if G were 0, then we would have 

Fo{x)F[{x)=F^{x)Fl,(x). (2.15) 

Now, let / C O be an open interval in O. If either F^^ or Fi were identically on /, our determinant would be 
identically on P (contradiction); it follows that we may choose / so that both and Fi are nonzero. Then 



we can divide both sides of ( 2.15 ) by Fq{x)Fi{x) and integrate; we find that Fq oc Fi on /. But this again 



(21 

makes the determinant 0. We conclude that the linear dependence satisfied by Vp must take the form 

V^p\yU = GV^^\y)2i (2.16) 

with G 7^ 0. 

f2') (2'] 

In particular, we find that the 2-form V' orthogonal to Vp (y) is not itself in the span of Vp (y)- In 
particular, any 2-form can be written as a linear combination of V and some of the Vf'{y). Taking the inner 

(2) 

product with Vp {x), we conclude that for < i < j < 3, we have 

V^'\x)., = p.,{x) (2.17) 

(2) 

for some polynomial p of degree at most 4. Now, since Vf{x) a Vf'{x) = 0, we find: 

Fo{x)pi2ix) - Fi{x)po2{x) + F2{x)poi{x) = (2.18) 
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for alH, j, fc. Similarly, since 



, we have 



Now, since V^\x)ki is linearly independent of Vp^'{x)ij, we can solve these two equations for Fi and as 
rational multiples of Fq, subsitute into the equation (a;) 12 = ^12(2;), then solve for Fq. We find: 

p P01P02 - P02P01 

-Tfl = 



^V^^\x) = Vf{x)AVI^{x) 



Fo{x)p[2{x) - Fi{x)p'o2{x) + F2{x)p'o^{x) = 

.(2) 



(2.19) 



(2.20) 



Fi = 



D 

Poip'i2 - Vi2Vm 



P02P'i2 - -P12P02 



where 



D = det 



Poi P12 P20 
P'oi P12 P20 



(2.21) 
(2.22) 

(2.23) 

(2.24) 



\Poi P12 P20 1 

Wc observe that each numerator has degree at most 6, as does the polynomial D. In particular, if we exclude 
any given _F, wc can express the squares of the other F as rational functions with common denominator of 
degree at most 6. It follows that the functions have at most 8 poles between them, and thus that we can 
write 



i^o(x) =Po(x)p(x)-l/2 

^^i(x) =pi(x)p(x)-i/2 

i^2(x) =P2(x)p(x)-l/2 

Fi{x)=Vi(x)v{x)-^^^, 



(2.25) 
(2.26) 
(2.27) 
(2.28) 



where po, Pi, P2, and are polynomials of degree at most 7 and p is a polynomial of degree at most 8. 

We now need to show that, in fact, each Fi is a polynomial of degree at most 3. By the usual factorization, 
we find: 



det(pj(a;) plAx) pj{y) p'Ay))o<j<i cc p{x)p{y){x - y)'^, 



(2.29) 



valid on M. Without loss of generality, we may assume that the constant of proportionality is 1, and that 
p{0) = 1. Dividing both sides by {x — y)'^ and taking the limit as a;, y — > 0, we find: 

/po(0) Pl(0) P2(0) P3(0)\ 
p'oiO) p'i(O) P'2(0) p'siO) 

p'm p'i'(o) p'm p'm 

\p'^'{0) p'{'{0) p'^'iO) p'i'{Q)j 



det 



(2.30) 
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Applying a suitable linear transformation to the polynomials pi, we have, without loss of generality, 



Poix) = 1 +po4:X^ +Pa5X^^ +Pa6x'^ + Pmx^ (2-31) 
Pi{x) — X + piix^ + Pibx''^ + piqx^ + pnx"^ (2.32) 

P2{x) = x'^ +P2ix'^ +P25X^ +P26x'^ + P27x'^ (2.33) 
P3{x) =X^ + P34X^ + P35X^ + P3fiX^ + psix'' (2.34) 

We can then solve for p{x) by taking y = above; we find 

p{x) = x~^{p2{x)p'^{x) - p3{x)p'2{x)). (2.35) 

At this point, we can compare coefficients on both sides of ( 2.29| ), obtaining a number of polynomial equations 
relating the coefhcients Pij , < i < 3, 4 < j < 7 . The resulting ideal can be verified (using magma, for instance) 
to contain the polynomials (p25 +P34P35 -Pse)^, (P26 +P34P36 -Pa?)^, and (p27 +^34^37)^; passing to the radical, 
we can then solve for pij^ < i < 2, 4 < j < 7. Substituting in, we find that p{x) is now a square, and that 
each Pi{x) is a multiple of y^p{x). In other words, each Fi{x) is a polynomial of degree at most 3, and we are 
done with the case 1 = 2. 

It remains only to show that we can reduce the cases Z > 2 to cases of lower dimension. Choose a particular 
element xq G O. By replacing the Fi with appropriate linear transformations, we may assume 

vP{xo)=CeiAe2, (2.36) 

for some nonzero constant C. In particular, we find that 

det{FJ{xi))o<^<l, o<j<2i = C'dct(Fj(xi))i<i<i, 2<j<2i (x G{xo)A{xi, X2, ■ ■ ■ xi^i)'^ Y[ xfGixi)"^. (2.37) 

l<i<l 

By induction, it follows that for 2 < i < Z, there exist polynomials Pi{x) of degree at most 21—1 such that 
Fi{x) = (x — xo)'^G{x)pi{x) on O. Undoing our linear transformations, we find that for every polynomial p{x) 
of degree at most 21 — 1 vanishing to second order at xq, we can write G{x)p{x) as a linear combination of the 
Fi{x). But this was independent of our choice of xq. In particular, taking x'q to be any other element of O, we 
have 



^ ^ 3(xo - x'q){x - xof + 2{x - x^f + 3(xo - x\){x - x'^f - 2{x - x'^;)- 

{xo - x'qY 



(2.38) 



and 



-2(xo - x'q){x - x^f -{x- Xpf - (xq - x'q){x - x'^f + {x- x'^f 
{x xo)- {xo-x',y • ^^-"^^^ 

It follows that for any polynomial p{x) of degree at most 21 — 1, G{x)p{x) is a linear combination of the Fi{x). 
By dimensionality, it follows that each Fi{x) is itself of the form G{x)p{x), and we are done. □ 
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3 Linear fractional transformations 



It will be convenient in the sequel to determine how matrix ensembles behave under a linear fractional change 
of variables. To be precise, let / be a weight function, and consider what the density of one of its associated 
matrix ensemble is in terms of the variables yi defined by Xi = {ayt + P)/{'yyi + 6). Clearly, we need only 
determine how A{x) and Hi dxi transform. 
We readily compute: 



aS — (3^ 



thus answering that question. As for A: 



Lemma 3.1. Let yo, yi, ... yi-i he a collection of I real numbers. Then for any a, (3, 7 and 6 such that 'yyi + S 
is never 0, 

^ f = - ^^y^'-"^'" W^y^ + ^f-'^ivi)- (3-2) 

\iyi + o J 

Proof. For each i < j, we have 

ayj + P ayt+ P a5 - (3^ 



M-Vi)- (3-3) 



Multiplying over i < j, we are done. □ 

We thus obtain the following transformation rules: 
Theorem 3.2. Let f be any weight function. Under the change of variables Xi = {ayi + P) / {-fyi + 6) , we have: 

OEiifix)) ^ 0E;(|a(5 - l3^\^'+'^/^{jy + 6)-'-'f{y)) (3.4) 
VEiifix)) ^ \JEi{\aS - H\^y + 5)-^'f{y)) (3.5) 

SEi{f{x)) ^ SEi{\a5 - h?'-\iy + 5f-^'hy)), (3.6) 

where f{y) = f{{ay + (3)/{'yy + S)), the normalization constants are the same on both sides, and for OE21, 
-yy + 6 must be positive over the support of f. 

Proof. When a6 — p'j < 0, the LFT reverses the order of integration, thus justifying the extra factor of (—1)' 
introduced for UE; and SE;. For OE;, there is a more subtle difficulty, namely that the relative order of the 
eigenvalues is significant, and can change. If we simply reverse the order, this is not a problem (the total effect 

is (— thus cancelling out the sign ofaS — (3j). So we can restrict to the case ad — > 0. The effect of 
the LFT is then to cyclically shift the ordering, taking the eigenvalues with x > a./"/ and making them smallest. 
If there arc k such eigenvalues, the sign of the Vandcrmondc matrix is changed by (— l)*^'-'"^'' ; thus if I is odd, 
there is no problem. On the other hand, if I is odd, we have a problem unless the eigenvalues are restricted to 
only one side of a/j, or equivalently that jy + 6 has constant sign over the support of (/). Since 

ay + P ^ -ay-P 
ly + o -7y - 6 

we may take this sign to be positive. □ 
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Remark. For algebraic purposes, we can often ignore the constraint "fy + S > 0, since the transform still has 
the correct form to be a matrix ensemble density, despite not being nonnegative. 

The upshot of this is that we can use this freedom to send a suitably chosen point to oo, thus simplifying 
our analysis below. 

4 The main results 

For a matrix ensemble Af , we define even(M) to be the ensemble obtained by taking the 2nd largest, 4th largest, 
etc. eigenvalues of M, and similarly for odd(M). 

When considering even(M) or odd(M) for M — 0E„ U 0E„ or AI = 0E„ U OE„-|_i, the following lemma is 
crucial: 

Lemma 4.1. For any integer n > 0, 

A{xs)A{x{o,i,...2n~i}~s) = 2"A(x{o^2,...2«-2})A(a;{i,3,...2„-i}) (4.1) 

SC{0,l,...2n-l} 

\S\—n 

and 

^ A{xs)A{x[o,i^,„2n}-s) = 2"A(a;{o,2,...2«})A(a;{i,3,...2„_i}) (4.2) 

SC{0,l,--2ra} 
|S|=n+l 

Proof. Consider what happens when we exchange Xi and XiJ^2 hi a term of either equation. If i, i + 2 G or 
i,i + 2 ^ S, then 

A(xs)A(x{o,i,...i-i}-5) ^ -A(a;5)A(a;{o,i,...i-i}-5), (4-3) 

since A is alternating. Otherwise, we see that every factor Xj — Xk with j > k is taken to another such factor, 
except for the factor Xi+i — Xi or Xi+2 — Xi+i, whichever is present. So each term in our sum is taken to the 
negative of a term from our sum; it follows that the sum is alternating under parity-preserving permutations. 
It follows that it must be a multiple of 

A(a:{o,2,...2r//2i-2})A(x{i^3,...2[;/2j-i})- (4.4) 

By degree considerations, it remains only to verify the constant, which we can do by considering the coefficient 
of largest degree in xq, and applying induction. □ 

Remark. The even case of this lemma is implicit in [|l^, where it was used to analyze even(OE„ U 0E„) with 
respect to the weight function 1 on the unit circle. 

From the lemma, it follows that the density of 0E„(/) U 0E„(/), expressed in terms of ordered variables, 
is proportional to 

W /(a;i)A(a;{o,2,...2n-2})A(a;{i,3_...2„-i}); (4.5) 

0<i<2n 
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similarly, the density of 0E„(/) U OE„_|_i(/) is proportional to 

Y[ /(2;i)A(x{o,2,...2n})A(a;{i,3,...2ri-l}). (4.6) 

0<j<2n 

For some weight functions /, if we integrate over the odd/even variables, the resulting density is the density of a 
unitary ensemble; we wish to determine precisely when that is. We first consider the case even(OE2(/)UOE3(/)). 

Theorem 4.2. Let f : R R be a function which is differ entiable on a possibly unbounded open interval 
/ C M and elsewhere. Suppose even(0E2(/) U 0E3(/)) = UE2((7) for some function g. Then up to a linear 
transformation of variables, f must have one of the following forms. On the interval (0, 1); 

f{x) (xx''{l-xf{l-rx)-'^-°'-'^, a,(3>-l,r<l (4.7) 

/(a;) oc 2;-^-"e-i/^(l -x)", a > -1. (4.8) 

On the interval (0, oo); 

f{x) cx x"(l - rxf, a>-l, a + /? < -3, r < (4.9) 

/(x) cx x-*^"e-i/^, a>-l, (4.10) 

/(x)cxx"e-^, a>-l. (4.11) 

Finally, on the entire real line: 

/(x) cx (1 + a;2)", a<-3/2 (4.12) 
fix) cx e-"'/2. (4.13) 

Proof. We need to integrate this over the variables X2i, and thus need to evaluate the determinant 



det(/ fix)x3dx)Q<,^j<2, (4.14) 

J [x2i-l,X2i + l] 

where we take x^i = a to be the left endpoint of /, and x^ — b to be the right endpoint of /. In particular, we 
need to determine when there exists a function g{x) with 

det( / f{x)x^ dx)o<.,j<2 oc g{xi)g{x3){x3 - xi). (4.15) 

J[a;2i-l,a;2i+i] 

As in section 10.6, we may use row operations to transform this to: 

det{FJ{x2^+l))o<^,J<2), (4.16) 

where we define 

Fjiy)= I f{x)x^dx. (4.17) 



We cannot quite apply theorem 2.1, however, since the last column of our determinant is constant. However, 
we clearly have -Fb(^) > Oj so we can eliminate that column, obtaining 



Fo{b) Aet{F,{x2^+i) - ^^Fo{x2^+i))o<.,.,<i). (4.18) 
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This, then, satisfies the hypotheses of theorem there thus exist linear polynomials pj such that 

P2{x){Fi{x) - CiFo(x)) = pi{x){F2{x) - C,Fo{x)), (4.19) 

where we have set 

Q = #!• (4-20) 
Fo{b) 

Differentiating twice and using the definition of Fi, we find: 

(P2(a:)(a: - Ci)~piix){x^ ~ C2))f'ix) = (~2(x - ^1)^2(2;) + 2{x^ - C2)p[ix) ~ P2ix) + 2xpi{x)) f (x) . (4.21) 

We can thus solve this for f'{x)/ f{x); we find that f'{x)/f{x) has the form p{x)/q{x) with deg(p) < 2, 
deg(q) < 3, and deg(a;p + 4g) < 2. We observe that these conditions are, naturally, preserved by linear 
fractional transformations. In particular, by applying a suitable linear fractional transformation, we may insist 
that q be strictly cubic, and that both endpoints of / be finite (possibly equal). (The result may very well no 
longer be a matrix ensemble, but as we noted above, this does not affect any algebraic conclusions.) 

Now, consider how f{x) and q{x) must behave at and 1 Differentiating ( [l.lS| ) once and taking a limit 
X X-i we find, since each Fi{x-i) — 0, 

lim ip2{x){x-Ci)+pi{x)iC2-x))fix) ^0 (4.22) 

X—>-X-i 

But this is just liuix^x-i Q{x)f{x). If q{x^i) ^ 0, then we must have liuix^x-i f{x) = 0. Then 

lim ^=00. (4.23) 

The only way this can happen is if q{x^i) = after all. Similarly, we have q{x2n+i) = 0. 

Suppose first that a ^ 6. Then up to LFT, we may insist that a — Q and 6=1, and thus q{Q) = q{\) = 0. 
We thus have two possibilities. The first is that q{x) has an additional zero, neither nor 1. In this case, 
integrating /'// and taking into account the constraints on p(x), we obtain 

J{x)=x°'{l-xf{l-rx)-^-"^f^ (4.24) 

Now, for / f{x) not to diverge at 0, we must have a > —1, and similarly P > —1. But then —4 — a — /3 < —2; 
it follows that (1 — rx) must be nonzero on (0, 1); in particular, r < 1. The other possibility is that q{x) has a 
double root, WLOG at 0. Upon integrating f / f, we obtain 

fix) = a;-"*-"e-''/^(l - x)", (4.25) 

and find a > —1, f3 > 0. The possibilities for (0, oo) then follow by LFT. 

The other case we must consider is / = R, and thus deg(g) = 2, deg(p) = 1. If g had a simple root in R, 
it would have two (WLOG and 1), and thus / would have the form a;"(l — x)^ with a,/3 > —1. But then 
the integral for F2 would diverge. Similarly, if q had a double root at 0, / would have the form x"e^^^^ , which 
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would diverge on one side of 0. Thus either q has a pair of complex roots, or deg(g) = 0. In the first case, a 
linear transformation takes the roots to ±i, and thus 



(4.26) 



For F2 to be; wcll-dcfincd, we must have a < —3/2. The other possiblity gives log(/(a;)) = ax^ +bx + c; thus a 
linear transformation gives 



fix) 



(4.27) 
□ 



We can now extend this to n > 2. 



Theorem 4.3. Fix an integer n > 2. Let / : M — > M fee a function which is differentiable on a possibly 
unbounded open interval 7 c R and elsewhere. Suppose even(OE„(/) U OE„+i(/)) = UE„(5f) for some 
function g. Then up to a linear transformation of variables, f and g can have precisely the following forms. On 
the interval (0, 1); 



f{x) (xx"{l- xf{l - ra;)-"-^-"-^, a, /3 > -1, r < 1 
g{x) oc x^'^+^l - xfl^+\l - rx)-2"-3-2«-2/3 

/(x) (xa;-"-2-"e-i/2^(l-a;)", a > -1 
gix) oc a;-2"-2-2«e-i/^(l - xf"+K 



On the interval (0, 00); 



f{x) oc x"{l - rxf, a> -1, a + (3 <-n-l, r < 
g{x)(xx'^"+^{l-rxf'^+^ 



fix) cx x-"-2-«e-i/2a.^ 
gix) oc a;-2"-2-2"e-i/^ 



a > -1 



fix) oc a;"e-^/2, 
gix) oc x^"+'e-' 

Finally, on the entire real line: 



a > 



/(x) oc (l + x^)", a<-(n + l)/2 
gix) oc (l + a;2)2«+i 



fix) 



gix) oc e ^ 



(4.28) 
(4.29) 

(4.30) 
(4.31) 

(4.32) 
(4.33) 

(4.34) 
(4.35) 

(4.36) 
(4.37) 

(4.38) 
(4.39) 

(4.40) 
(4.41) 
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Proof. As for n = 2, the issue is when 

Aei{Fj{x2i+i))o<i,j<n) (4.42) 

takes the form of an orthogonal ensemble. Applying an LFT as necessary, we may assume that a = —oo. Now 
differentiate with respect to xi, divide by x"~^f{xi), and take a Umit as xi — > — oo. On the one hand, this 
operation takes orthogonal ensembles to orthogonal ensembles. On the other hand, we can then expand along 
the first column, finding that 

det(F,_i(x2i+i)) ) (4.43) 

must take the form of an orthogonal ensemble. By induction, we find that f{x) must satisfy the constraints 
valid for n—1. Upon undoing the LFT, we obtain the desired "only if" result. 

It remains to show that each of the above weight functions actually do work. We need only consider the 
following possibilities: 

f{x) =x°'{l-xf (4.44) 
f{x) = a;"e-^ (4.45) 

f{x) = {l+x^r (4.46) 
f{x) = e-'l\ (4.47) 

(on (0, 1), (0, oc), K, and K respectively) since the others are all images of these under LFTs. 

For j{x) = a;"(l — .x)''*, observe that 

(q + /3 + j + 2)F,+i(a;) - (a + j + \)Fj{x) = -x\x'^^^{\ - xf^^); (4.48) 

this is true for a; = 0, and both sides have the same derivative. In particular, for each j, we have a polynomial 
(a;) of degree i — 1 and a constant Cj with 

F^{x) = Qi^o(a;) +p,(a;)(a;«+^(l - xf^\ (4.49) 

In particular, this must be true for x = 1, and thus Cj = Fj{l)/Fo{l) as required. We thus find that we obtain 
a unitary ensemble with weight function proportional to a;^"+^(l — x)'^^~^^. Similarly, for f{x) = a;"e~^, we 
have 

Fj+i {x) - (a + i + l)Fj {x) = -a;«+^+ie-^, (4.50) 

so g{x) oc x^"+^e~^^. 

For f{x) = (1 + a;2)«, we find 

(2a + j + 2)F,+i(x)+jF,_i =xJ'(1 + .t2)"+i (4.51) 

This allows us to solve for each Fj except Fq; we obtain g(x) oc (1 + x^)^""*"^. Finally, for f{x) = e~^^/^, we 
have: 

F,•+l-jF,•_l=-x^e--'/^ (4.52) 

and ^(a;) a e~^^ . □ 
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Remark. We observe that in each case g{x) oc f(x)'^q{x). 

For even(OE„ U 0E„), the calculations are analogous, and we have: 

Theorem 4.4. Fix an integer n > 2. Let / : R — s- M fee a function which is differentiable on a possibly 
unbounded open interval / C M and elsewhere. Suppose even(OE„(/) U 0E„(/)) = UE„((7) for some function 
g. Then up to an order-preserving linear transformation of variables, f and g must have one of the following 
forms. On an interval with right endpoint 0; 

/(x) cx (-a;)"(l -rx)-"-!-", a > -1, 1 ^ r/ (4.53) 
gix) cx - ra;)"2"-i-2" (4.54) 

fix) cx (-x)-"-iei/^ (4.55) 
gix) oc (-.x)-2»e2/". (4.56) 

On an interval of the form (a, 00), a > —00: 

/(x) cx (1-ra;)", a < -n, r < (4.57) 
g(a;) cx (1 - ra;)2"+i (4.58) 

fix) cx e-^ (4.59) 
gix) cx (4.60) 

On the entire real line, no possibilities exist. 

Remark. The relation between g and / is here slightly modified, by removing the factor of q corresponding to 
the left endpoint; similarly, for odd(OE„ U 0E„), we remove the factor corresponding to the right endpoint, 
and for odd(OE„_i U 0E„), we remove both factors. 

For odd(OE„ U 0E„), we need simply reverse the ordering. For odd(OE„ U OEn+i), we have: 

Theorem 4.5. Fix an integer n > 2. Let / : M — > M fee a function which is differentiable on a possibly 
unbounded open interval / C M and elsewhere. Suppose odd(OE„_i(/) U 0E„(/)) = UE„(5) for some 
function g. Then f and g have the form 

/(x) cx (1-ra;)-", (4.61) 
g(a;) cx (1 - ra;)-2"+i (4.62) 

for some r (possibily 00) with l/r ^ L. 

Proof. The only tricky aspect of this case is that the determinant we must analyze is no longer of the form to 
which Theorem 2.1 applies; to be precise, we need 

detiFjix2i+2) - Fjix2i))a<i<„ (4.63) 
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to have orthogonal ensemble form. But this determinant is clearly equal to the determinant of the block matrix 

( ^'^^<-^<- ] (4.64) 



Adding the first column to the other columns, we can then apply Theorem 2.1, and argue as above. □ 

We finally consider a fifth possibility for decimation. Recall that for the circular ensemble results cited 
above, while there was a local notion of order, there was no notion of largest. This suggests that we consider the 
ensemble derived by choosing randomly between odd(M) and even(M). More precisely, for an ensemble with 
an even number of variables, we define alt(M) to be even(M) with probability 1/2 and odd(M) with probability 
1/2. 

Theorem 4.6. Fix an integer n > 2. Let / : M — s- R 6e a function which is differentiable on a possibly 
unbounded open interval / C K and elsewhere. Suppose alt(OE„(/) U 0E„(/)) = UE„((/) for some function 
g. Then up to a linear transformation of variables, f and g have the form 

/ = (1 + .t2)-("+i)/2 (4.65) 
<7 = (l + x2)-". (4.66) 

Proof. Consider the determinants associated to even(OE„(/) U 0E„(/)) and odd(OE„(/) U 0E„(/)). Up to 
cyclic shift, only one column differs between the two determinants, thus allowing us to express their sum as a 
determinant. When n is even, the 'special' column takes the form 

(F.ixo) + F,{x2n-l) - Fj{I))o<,<n; (4.67) 

here Fj{I) ~ J^^jX-' f{x). Taking appropriate linear combinations, we obtain the determinant 

det(i^,(x,) - i^;(/)/2)o<^j<n. (4.68) 
When n is odd, the special column takes the form 

iF,{x2n-l) - Fjixo) - Fj{I))o<,<n; (4.69) 

this leads (up to sign) to the n + 1 x n + 1 block determinant 

detf° (^^(^))"^^<" V (4.70) 

\1 {F.j(xi))o<ij^n) 

We first analyze the case n odd. In this case, the usual theory tells us that there exist polynomials Pj{x) 
and q{x) of degree at most n — 1 with 

F,{x)-C,F^{x)^p,{x)/q{x) (4.71) 



17 



for all j, with Cj as above. Now, evaluating this at an endpoint of /, wc find that the polynomials pj must 



have a common root (possibly oo). In particular, it follows that / must satisfy the conditions of Theorem L3 
On the other hand, we find that 



dx q{x) 



x'f{x)-C,J{x)^^^ (4.72) 



for each j; in particular, f{x) must be a rational function. We therefore have the following possibilities to 
consider: 

f{x) = x"(l - xf, a, /? e N, (4.73) 

/(x) = (1 + a;2)"", aeN,a>n/2, (4.74) 

on (0, 1) and M respectively. In the first case, we find that each 

F,{x)-C,Fq{x) (4.75) 

is a polynomial of degree j + 2. In particular, Fn-i{x) — Cn-iFQ{x) is a polynomial of degree contradicting 
the bound on Aeg{j)j[x)). 

In the second case, we observe that 

F,{x) - C,Fo{x) = r,ix){l + x')'-" (4.76) 

for polynomials rj{x) of degree j — 1. In particular, we find that pi{x) / q{x) cx (1 + x^)^^", implying, since 
a > n/2 > 1, that 

g(x) cx (1 +.t2)"-i (4.77) 

Since deg(q) < n — 1, we have n/2 < a < {n + l)/2, the only integral solution of which is a ~ (n + l)/2. In 
this case, the relevant degree bounds all hold, and thus the determinant is indeed of the correct form, giving 
g{x) oc (1 + x^y^'^" = (1 + x^)~^ as required. 

For n even, we must have polynomials pi of degree at most n — 1 with 

p^{xKF,{x) - C;) - Pj{xKF,{x) - CD, (4.78) 

where we write C'j for Fj{I)/2. We can rewrite this as: 

Po{x){F,{x) - C.Foixj) = {p,{x) ~ C,po{^)){Fo{x) - C^), (4.79) 

using the fact that Cj = C^/Cq. For n > 2, we conclude that f{x) must satisfy the conditions of Theorem 
4.3. Now, if the endpoints of / are different, then we find, since Fo{x) — Cg = ±Cq at both endpoints, that 
each Pj{x) — CjPo{x) = at both endpoints. But this causes the polynomials to be linearly dependent, a 
contradiction. On the other hand, in the other cases, we know that C'l — and pi cx 1. In both cases, we obtain 
from the identity for Fi{x) — CiFo{x) a differential equation for pq{x). For e^^ no polynomial solution to 
the equation exists. For (1 + a;^)~", we can find an explicit power series solution to the equation, and find that 



18 



a polynomial solution exists only when a is half-integral, when the solution has degree 2a — 2. As above, this 
leaves only one possibility for a, namely a = {n+ l)/2, as required. 

It remains to consider n = 2. Here we can twice differentiate the equation 

Po{x){F,{x) - C[) = pi{x){Fo{x) - C^) (4.80) 

(with po and pi linear) to deduce that 

f{x)/f{x)=pix)/q{x) (4.81) 

with deg(p) < 1, deg{q) < 2, and deg{xp+ 3q) < 1. So up to LFT, f{x) must have one of the forms 

f{x) =x", a> -3/2 

fix) = e-- 

f{x) = {l+x-')-'/\ 



(note that if we exchange and oo in the first case, we replace a with —3 — a, justifying our restriction on a.) 
In the first case, if a ^ — 1, the following: 

(4-86) 

must be a linear polynomial for suitable constants Di and Dq respectively proportional to C[ and Cq (and thus 
Do ^0). We deduce therefore that a = 0. But then we readily determine that only the empty interval satisfies 
the requirements. For x~^ and e~^, there are not even appropriate choices for Do and Di (since both log(a;) 
and are transcendental functions). Finally, for the third choice, we readily verify that decimation indeed 
works as required. □ 

We now turn our attention to decimations of single orthogonal ensembles. We have, quite simply: 

Theorem 4.7. Fix an integer n>2. For any functions f and g with f dijferentiable on a possibly unbounded 
open interval 7 c M and elsewhere, each of the following pairs of statements is equivalent: 



cvcn(OE2„(/)UOE2„+i(/)) 


= UE2„(5) 


and 


even(OE2„+i(/)) = SE„((5//)2) 


(4.87) 


even(OE2„(/)UOE2„(/)) 


-UE2„(.g) 


and 


even(0E2„(/)) = SE„((.g//)2) 


(4.88) 


odd(OE2„(/)UOE2„(/)) 


-UE2„(5) 


and 


odd(OE2„(/)) = SE„((5//)2) 


(4.89) 


odd(OE2„-i(/)UOE2„(/)) 


= UE2„(5) 


and 


odd(OE2„_i(/)) = SE„((g//)2) 


(4.90) 


alt(OE2„(/)UOE2„(/)) 


= UE2„(5) 


and 


alt(0E2„(/)) = SE„((5//)2). 


(4.91) 



Proof. Consider, for instance, even(0E2„+i(/)). Once we integrate along the largest, 3rd largest, etc. variables 
and do some simplification, the resulting matrix has columns {Fj{x2i+i))o<j<2n and {x2i_^_if{x2i+i))o<j<2m 
with the last column given by {Fj{b))o<j<2n- In particular, we note that aside from the last, constant, column. 
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the columns come in pairs, one the derivative of the other. Thus the determinant is essentiaUy of the form 
considered in Theorem 2^ . In particular, it has a symplectic ensemble form if and only if the determinant 



det(Fj(a;,+i))o<ij<2« (4.92) 

(of which our determinant is a derivative) has an orthogonal ensemble form. But by the proof of Theorems 
4.2 and 4.3, this precisely what we needed to show. (The statement about the resulting weight fimctions is 
straightforward.) Similar arguments apply for the remaining equivalences. □ 



5 Random matrix applications 

In random matrix applications /, g and {g/ f)"^ must be (up to the scale of x) one of the four classical forms 
(1.12) and ( 1.16| ). So specializing Theorems 47^4.5 we can read off for which of the classical forms the statement 
of Theorem 4.7 is valid. 

Theorem 5.1. Restricting attention to the classical weights (l.li) and the statement (Jf-Sl) holds for 



(/,5) = 



(e-^'/2^e-^') 

(2,(a-l)/2g-x/2^^ag-:.)^ a; > 

(2.(a-i)/2^ _ a;)(''-i)/2,a;'^(l - xf), < X < 1 
((l + a;2)-("+i)/2,(l + x2)-"). 



the statement {4 -St) holds for 



if, 9) 



(e-^/2,e-^) 

((1 - x)('^-i)/2^ (1 - x)""), < x < 1 



while (4-.8!> ) is valid for the particular pair of Jacobi weights 



(5.1) 



(5.2) 



(5.3) 



and (4-9C) is valid for the particular pair of Jacobi weights 

(/,5)-(l,l). 



(5.4) 



Because the weights in Theorem 5.1 occur in the matrix ensembles listed in the Introduction, the theorems 
of Section imply inter-relationships between the different ensembles. 



Theorem 5.2. The following relations hold between the above matrix ensembles under decimation, for all 
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n > 0; 



even(G0E2„+i) 


— GSE„ 


(5.5) 


even(Synim(n; C)) 


= Anti(n; C) 


(5.6) 


even(Mat(2p+ \ ,2q+\-K)) 


Mat(p, g;H) 


(5.7) 


even(Beta(2pi + 1, 2p2 + 1, 2g + 1; M)) 


= Beta(pi,p2,g;]3) 


(5.8) 


eveii(GOE„UGOE„+i) 


= GUE„ 


(5.9) 


even(Symm(n; C) U Symm(n; C)) 


= Mat(n,n;C) 


(5.10) 


even(Symm(7i; C) U Symm(rt + 1; C)) 


= Mat(n + l,?i;C) 


(5.11) 


even(Mat(p,9;R) UMat(p+ 1;M)) 


= Mat(p, g;C) 


(5.12) 


even(Beta(pi,p2,9;IR) U Beta(pi + l,p2 + l,q + 


= Beta(pi,p2,g;C) 


(5.13) 



Remark 1. It would be very nice to have a direct, matrix-theoretic, proof of any of the above relations. 
Remark 2. There are actually a few more relations, all of which follow from the above together with the relation 

Mat(n + l,n;R) = Symm(n;C). (5.14) 
Again, a matrix-theoretic proof of this would be nice. 



We now turn our attention to the implications of Theorem 5.1 with respect to gap probabilities. In circular 
ensemble theory the results (1-14) and ( 1.15| ) were shown |^ to imply inter- relationships between the probability 
of an eigenvalue free region amongst the various symmetry classes. With E''^'^ {p; J; n) denoting the probability 
that, for the ensembles COE„ {(3 — 1), CUE„ {f3 — 2) and CSE„ (/3 = 4), there are exactly p eigenvalues in the 
interval J, the inter-relationships are 

(0; (0, s);n) = E^'^ (0; (0, s); n) [e^'^ (0; (0, s);n) + E^''^ (1; (0, s);n)) 

E^^\p; (0, 5); n) = E^'\2p; (0, s); 2n) + ^E^'H^P - 1; (0, s);2n) + ^E'^'\2p + 1; (0, s);2n) (5.15) 

where 

E'-'^\p]J;n):^0, for p<0. (5.16) 

Similar inter-relationships between gap probabilities, but now with the eigenvalue free interval including an 



endpoint of the support of the interval, can be deduced from the pairs of statements of Theorem 4.7 



Theorem 5.3. Let E^^\p\ J\ g\n) denote the probability that, for the ensembles 0E„(5) (j3 — 1), \JEn{g) 
((3^2) and SEn{g) ((3 — 4,) the interval J contains exactly n eigenvalues. The statements { 4-8'\ ) imply 



i;(2)(0; J-g- 2n) = £;(i)(0; J; /; 2n)£;(i)(0; J; /; 2n + 1) + £;(i)(0; J; /; 2n)£;(i)(l; J; /; 2n + 1) 
W^HO; J; /; 2n -f l)£;(i) (1; J; /; 2n) 
E'^^\p;J-{g/ff-n) = ii;(i'(2p;J;/;2n+l)+i?W(2p-f-l;J;/;2n-t-l) (5.17) 
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for J — (—00, — s) or (s, 00); the statements {4-.8t ) imply 



i?(2)(0; -s);g; 2n) = (i?(i)(0; (s, 00); /; 2n))' + 2i?(i)(0; (s, ^); /; 2n)i?(i)(l; (s, /; 2n) 

£;(4) (p; (5, 00); n) = E^^^ (2p; (s, 00); /; 2n) + (2p + 1; (s, 00); /; 2n) 



(5.18) 



and 



i?(2)(0;(-oo,-5);g;2n) = 
i?(4)(p;(-cx3,~s);(5//)2;„) = 

the statements ( 4-8!\ ) imply 

i?(2)(0;(s,^);5;2n) 



(i?W(0;(-(X3,-5);/;2n)^ 
i?(i)(2p; -s); /; 2n) + iJ^) {2p ~ 1; (-ex., -s); /; 2n); 



(5.19) 



i?«(0;(s,oo);/;2n)^ 

(2p; (s, 00); /; 2n) + E<^'^ (2p - 1; (s, ex.); /; 2n) 



(5.20) 



and 



£;(2)(0;(-(X3,-s);5;2n) 

= (0; (-00, -s); /; 2n))' + 2^(1) (0; (-00, -s): /; 2n)^(i)(l; (-00, -s); /; 2n) 

i5(4) (p. -s); (5//)2; n) = E^'^ (2p; (-00, -s); /; 2n) + E^'^^ {2p + 1; (-(X3, -s); /; 2n); (5.21) 



f/ie statements (4-9C) imply 



S(2)(0; J;5;2n) = (0; J; /; 2n)s(i) (0; J; /; 2n - 1) 

- i?«(2p; J;/;2n- l) + i?«(2p- l; J;/;2n- l) 



(5.22) 



/or J = (— oo, —s) or (s, oo), while the statements ( 4. 91 ) imply the relations ( 5.1^ ) with n replaced by 2n in the 
first equation and (0, s) replaced throughout by J, J — (— oo, — s) or (s,oo). 



Proof. We will consider only the deductions from ( |4.87| ), as the other cases are similar. Let J be a single 
interval which includes an endpoint of the support of / and g. From the first statement in (4.87) we see that 
the event of a sequence of eigenvalues from \JF,2nig) not being contained in J occurs in three ways relative 
to the ensemble 0E2„(/) U 0E2„+i(/): (i) the eigenvalues from 0E2„(/) and those from 0E2„+i(/) are not 
contained in J; or (ii) one eigenvalue from 0E2„+i(/) is contained in J and no eigenvalue from 0E2„(/) is 
contained in J (note that the one eigenvalue must be either the largest (smallest) eigenvalue when J contains 
the right (left) hand end point); or (iii) one eigenvalue from 0E2,i(/) is contained in J and no eigenvalue from 
0E2„+i(/) is contained in J. This gives the first equation in ( 5.17 ). From the second statement in ( 1.87 ) we 
see that the event of a sequence of eigenvalues from SE„((f;//)^) containing p eigenvalues in J can occur in 
two ways relative to 0E2„+i(/): (i) there are 2p eigenvalues from 0E2„+i(/) in J; or (ii) there are 2p + 1 
eigenvalues from 0E2„+i(/) in J (ofwhichp+1 are integrated out in forming even(0E2ri+i(/))). This implies 
the second equation in ( 5.17 ). □ 
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Recalling (5.16) we see that in the case p = the equations (5.19), ( 5.20D and (5.22) give particularly 
simple inter-relationships between the i?('^)(0; . . . ). In fact referring back to Theorem ^.l| for the permissable 
pairs (/, 5) in these cases it is a simple exercise in changing variables to compute the E^^^O; . . . ) in terms of 
elementary functions. Recalling 



E'^^\0;J;w;n) := 



•' 1=0 0<j<k<n-l 



\Xk - Xjf, 



where J = (—00, 00) — J and C is such that E^^'' (0; 0; w;n) = l we find 
i;(i)(0;(0,s);e-"/2;n) = e-^"/2 



£;(2)(0;(0,s);e-^;n 
£;(i)(0;(0,s);(l-x)('^-i)/2;n 
£;(2)(0; {0,s);{l-xr;n 
E(^^\0;{0,s);{l~xr+';n 



= i;(^)(0;(0,s);e-";n) = e-"" 

= E''^\0; (1 - s, 1); n) = (1 - s)"("+'^)/2 

= i;(2)(0; (1 - s, 1); x"^; n) = (1 - s)"('^+«) 

= i;('')(0;(l-s,l);a;'^+^n) = (1 _ s)2«'+«« 



(in the first two cases the weight functions are restricted to x > 0, while in the remaining cases < a; < 1). 
The equations ( 5.19 ), ( 5.2C ) and ( 5.22 ) for p = can be checked immediately. 

The pairs of equations ( ^.17 ), (5.18) and (5.21) contain E^^\l] . . . ) as well as £^^''^(0; . . . ). In the equations 
( ^■18D and ( |5.2l| ) the dependence on i;(i)(l; . . . ) can be eliminated. Noting from Theorem |5.l| the allowed pairs 
(/, (7) for the validity of equations (5.18) and (5.21) the following result is obtained. 

Proposition 5.4. For (/, g) = (e'^/^, e""^), {x>0), and J = {s, 00), or for (/, g) = ((1 - x)^"-^'^/'^, (1 - x)"), 
(0 < a: < 1), and J = (1 — s, 1) we have 



.■'Ho;.;(.,//.V.,4(.''Ho;.;/;2„,.|S|^). 



(5.23) 



)(0;J;/;2n). 

In the scaled n 00 limit, as appropriate for the particular choice of weight function in (^.l[), the pair of 



equations ( 4.87 ) also imply an equation of the form ( 5.23 ). First consider the Gaussian ensembles with the 



scaling |11 



which corresponds to studying the distribution of the eigenvalues at the (soft) edge of the leading order support 
of the spectrum. Defining 

Eill{p;is,oo)) := l^^E^') (^p,i{2ny/^ + ^^j^,ooy,e-^'/'-,n) 
Eili{p;{s,^)) := lun^E(^^{p,ii2ny/^ + ^^,^y,e-^';n) 
Eiti{p;is,oo)) := ^lim i?(4)(p,((2n)V2 + _^i_^,^);e--V/2), 



(the existence of these limits is known from explicit calculation |22l; see below). The equations (|5.17| ) imply: 
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Proposition 5.5. For the scaled infinite Gaussian ensembles at the soft edge 

E^ms.oo)) = i(i.i^)(o; (.,00))+ (5.24) 

Eillil;{s,^)) = ii;it(0;(0,oo))-£;W(0;(s,oo)) (5.25) 

As alluded to above, the i?gQf|(0; (s,oo)) are known exactly from the work of Tracy and Widom To 
present these results, let q{s) denote the solution of the particular Painleve II equation 

q" = sq + 2g3 

which satisfies the boundary condition q{s) ~ Ai(s) as s cx3. Then we have 

i?l'f\(0;(s,oo)) = exp(^- J^^{t^s)q'{t)dt) 
(£;W(0;(s,cx3))j = E^ll{0;is,^))exp[- J q{t)dt) 

(£;(„t(0;(s,oo))) = EilliO;is,^))cosh^(^-J q{t)dt), (5.26) 

(in |2^ -E'ioft defined with s i— > s/2^/^ relative to our definition). The equation (5.24) is immediately seen to 
be satisfied, while the second equation gives 

(l;(s,cx3)))' =£;i'^,(0;(s,(X3))sinh2 (lj^q(t)dt). (5.27) 

Next consider the scaled limit at an edge for which the weight function is strictly zero on one side. For the 
classical ensembles this occurs in the Laguerre and Jacobi case; for definiteness consider the Laguerre case. The 
appropriate scaling is Q 

X 

and we define 

4ardfe(0>s);(«-l)/2) lim £;(i)(p;(0,s/4n);:r('^-i)/2e-^;n) 

n — >oo 

4a!-d(P; (0, s)-a) := lim E'^^\p; (0, 5/4n); x-^e"^; n) 

n — 'oo 

4ardfe(0. := lim £;(4)(p;(0,s/4n);a:'^+ie--;n/2) 

n — >oc 

(the existence of these limits for general a > — 1 can be deduced from the existence of the fc-point distributions 
in the same scaled limits [^). Use of (4.87) then gives the analogue of Proposition |5.5| for the hard edge. 

Proposition 5.6. For the scaled infinite Laguerre ensembles at the hard edge 

41(0; (0,.); a + 1) = ^(41(0; (0, {a - i)/2) + ..^V^' ^"'Z^' "\ ) (5.28) 
4ard(l;(0,s);(a-l)/2) = 4tid(0;(0,s);a+l)-42-d(0;(0,s);(a-l)/2). (5.29) 
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Exact Pfaffian formulas are known for eI^^^^{0; (0, s); (a — l)/2) and -Ejj^|.^(0; (0, s);a + 1) in the case a an 
odd positive integer while £'hard('^' (0, s);a) can then be expressed as a determinant (the dimension 



of the Pfaffians and the determinants are proportional to a), although ( 5.28| ) is not a natural consequence of 



these formulas. There are also multiple integral expressions for the same expression but again they do not 



naturally satisfy ( 5.28 ) 



6 Distribution functions for superimposed spectra 

In general, for a symmetric PDF p{xo, . . . ,a;„_i) the fc-point distribution function pk is defined by 



Pk{xo, . . . ,Xk-i) n{n- 1) ■ ■ ■ {n- k + 1) / p{xo, ■ ■ ■ , Xn~i)dxk ■ ■ ■ dxn-i- (6.1) 

In this section we take up the task of computing pk for even(Af ), odd(M), alt(M) with M — 0E„(/) U 0E„(/) 
and even(M), odd(M) with M = 0E„(/) U OE„+i(/). 
For M = 0E„(/) U 0E„(/) write 

D'=™"W(xo,...,x„_i) := det{F,{x,)-F,{I))o<^,J<n 
D°'"'('^\xo,...,Xn-i) ■■= det(F,(a;0)o<»j<n 
D""(^'(a;o,...,a;„_i) := det{Fj{x,) ~ ^FJ{I))o<^J<n (n even) 

(^;(/))o<,<« 

(l)l<i<ri-l {Fj{xi))o<ij^n 



nalt(M)/ X 1 ; / (-Pj(-^))o<J<« \ , 

D ^ '{xo, . . . ,Xn-i) := dct (nodd) 



and for M = 0E„(/) U OE„+i(/) let 



£»^™"(^)(a;o,...,x„_i) := det ( o<7<n+i 

{Fj{I))o<'j <n+l 
,.UT^,j\ I (l)o<j<n+l , 

i.-(-)(.o,...,.„_0 :^ det( (^^.(,^)) I (6.2) 

0<j<n+l 



In each case, workings contained in the proofs of Theorems 4.2 and 4.6 show (after relabelling the coordinates) 
that the PDF is proportional to 



Y[ f(Xi)A{xo, . . . , Xn-l)D{xo, . . . , Xn-l) (6.3) 

i=0 

for D as specified. Now introduce a set of functions {?7j(a;)}o<j"<n such that 

D{xo, . . .,Xn~i) oc det(77j(a;i))o<ij<„ 
and a set of monic polynomials {'7j(a;)}o<j<ji, Qj of degree j, such that the biorthogonality property 

f{x)qi{x)r]j{x) dx = Si^j 
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holds (assuming such biorthogonal famihcs exist). The fc-point distribution can be expressed in terms of these 
functions. 



(6.4) 



Lemma 6.1. For the PDF ^6.^\ ) and {ilj ix)}o< j<n, {Qj{x)}o<j<n specified as above, we have 



k-l 



Pkixo, . . . ,Xfc_i) = Yl fi^i) det ^ qi{xi)r]i{xj) 



0<i,j<k 



j=0 1=0 

Proof. From the definitions of {f7j(-'J;)}o<j<n and {qj{x)}o<j<n we see that ( [6.3^ is proportional to 

n-l 



Y[ /(a;i)det((7j(a;i))o<ij<„det(77j(xi))o<ij<„. 



The biorthogonal property allows the integrations required by the definition (OA) to be computed to give 
(U). □ 

Remark. Suppose for some {^j(a;)}j=o,....n-i we can write 

D{xo, . . . ,a;„_i) cx det(fj(a;i))o<jj<„-i. 
It is easy to show ^] that sufficient conditions for the existence of the biorthogonal sets is that 



det 



f{x)x'^j{x) dx 



0<i,j<p 



^0 



for p — 0, ... ,11 — I. 



For / a classical weight and so of the form ( 1.12 ) or ( 1.16 ), the biorthogonal functions can be computed ex- 
plicitly. This is possible because of the following special property of the classical weights and their corresponding 
orthogonality 



Lemma 6.2. Consider the pairs {f,g) of classical weight functions (5.1). Let {pj{x)}j^o,i,--- '^f 
monic orthogonal polynomials, Pj{x) of degree j corresponding to the weight function g, let {pk,Pk)2 denote 
their normalization with respect to integration over the measure g(x)dx, and define jk so that 



lk{Pk-,Pk)-2 



= < 



in the four cases respectively. With 



1 d g{x) 
f{x) dx f{x) ■ 



i(2fc + 2 + a + 6) 
a — fc — 1 



Ck -.^ Jk{Pk,Pk)2{Pk+l,Pk+l)2 



npk{x) = — —pk+i{x) + ^^^^ — —pk^i{x). 

[Pk+l,Pk+l)2 (Pk-l,Pk-l)2 
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Proof. This is a simple consequence of the property 



□ 



As a consequence, the determinant formulas (3.2) and (6.2) for D{xq, . . . , Xn-i) can be simplified. Let 



[u{x)]j := 



1=3 



{Pl,Pl)2 



J-1 



■pi{x) ^u{x)-Y^ 



and define 



rf\x) = 



1=0 



rf\x) = 



{u{x),pi{x))2 
{Pl,Pl)2 



Pi{x) 



(6.5) 



\ (4) / N ('■)-l.'---l)2 (2) . ^ 



Then by adding appropriate linear combinations of the columns in the determinant formulas (3.2) and (6.2), 
and making use of Lemma 3.2 we readily find for AI — 0E„(/) U 0E„(/) 





■ 1 ^n- 




(p-,(Xi)) 0<i<n 
0<j<n-l 




D°'''"^^'\xo,.. 


■ ■) •^n- 


-''=<n?;::5-( 


iPj{Xi)) 0<i<n 
0<j<n-l 




D-''(^'\xo,.. 


■ 1 -^n- 




{P]{Xi)) Q<i<n 
0<j<n-l 




D"i*(^^)(xo,.. 


■ 1 '^n- 




{Pj{Xi)) Q<i<n 
0<j<n-2 





(6.6) 



while for M = 0E„(/) U OE„+i(/) 



n-1 



1 -^n— 1 



g(a^») 

^ /(a;.) 



det 



0<'i,j<n 



TT ^detf (P,(a;0) o<.<„ {ri%{x,))o<^<n {rl^ll)o<^<n ) . (6.7) 

/(a^i) ^ 0<]<n-2 J 

In each case, setting ?7j(a;)/ f{x)pi{x)rjj{x) dx equal to g{x)/f{x) times the function in column j and 
^^(a;) = Pi{x) we have that 



f{x)qi{x)r]j{x) dx = Sij (i, j = 0, . . . , n - 1), 



which is the desired biorthogonality property. Hence substitution of these values into (6.4) gives the fc-point 
distribution in each case. 
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In particular with M = 0E„(/) U OE„+i(/) and / one of the classical weights in (p^) we read off that 



n-l 



This is the well known expression for in UE((7), and thus is in keeping with the result of Theorem 4.3, giving 
even(OE„(/) U OE„+i(/)) — UE„(g) for each of the pairs (/, g) in (5T). Furthermore, the Christoffel-Darboux 
formula evaluates the sum as 

S,{x,y) := y MeMM ^ (g(^).9(y))^/^P»(^K-i(^)-P«-i(^K(2/) (g g) 

{Pl,Pl)2 [Pn-l,Pn-l)2 X-y 

7 Distribution functions for alternate eigenvalues in a single 0E„ 

The /c-point distribution function for the alternate eigenvalues in a single 0E„ has a different structure to pk 
for the superimposed 0E„ spectra. The cases n even and n odd must be treated separately. 



n even 



Consider first even(OE„(/)) with n even. From the manipulations sketched in the proof of Theorem 4.7 we 
have that the PDF of this ensemble is given by 



- W f{x2i)Aet 
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1=0 



(7.1) 



0<i<n 
0<j<2Ti 



To perform the integration required by (BJ^) we introduce the skew inner product 



{u\v)i 



dxf{x)\u{x)l dyf{y)v{y)-v{x) I dy f{y)u{y) 

) ^ J X 

nOO 

dxf{x)u{x) / dyv{y)sgn{y~x), 



(7.2) 



together with a corresponding family of monic skew orthogonal polynomials {-Ri(a;)}i=o,i,... which are defined 
so that 



(i?2'i|^2j + l)l — —{Ii2j+l\R2i}l = fjSij, {R2i\R2j)l = (^2i+l|-R2j + l)l — 0. 

Note that the skew orthogonality property still holds if we make the replacement 

R2i+l{x) ^ R2i+l{x) +j2iR2tix) 



(7.3) 



(7.4) 



for arbitrary "f2i. However a Gram-Schmidt type construction shows {Ri{x)}i^o,i,... is unique up to this trans- 
formation. 



We will first express (7.1) as a quaternion determinant involving {i?i(2:)}i=o,i,... and then show how the 
property (7^) can be used to perform the integrations. This requires the definition of a quaternion determinant. 
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We regard a quaternion as a 2 x 2 matrix, and a quaternion matrix as a matrix with quaternion elements. With 
n even and 

Zn ■■= l«/2 I ^ Q I ' 

a n/2 X n/2 quaternion matrix Q is said to be self dual if 

In terms of its 2 x 2 sub-blocks this means that the quaternion element in position (kj) is related to the element 
in position (jfc), j < fc by 

f d -b\ f a b\ 

''^=[^c a ) '''' '^'^[c d)- 
Now for a self dual quaternion matrix the determinant, to be denoted qdet, is defined by |^ 

qdetO= {-^r^^''l[{Qabqbc---qdaf^ (7.5) 

where the superscript (0) denotes the operation ^Tr, P is any permutation of the indicies (1, . . . , n/2) consisting 
of / exclusive cycles of the form {a ^ b ^ c ^ ■ ■ ■ d a) and (— l)"/^-' is the parity of P. Furthermore, qdet Q 
is related to the Pfafhan via the formula ||] 

qdetQ = Pf QZ-\ 

which since {Pf QZ^^)^ — detQ (where here Q is regarded as an ordinary nx n matrix) implies [ p^ 

detQ = qdet(QQ^) (7.6) 

assuming det Q is positive. 

Proposition 7.1. With p{xo, X2, ■ ■ ■ ,Sn-2) denoting the PDF ([7.4), {i?i(x)}i=o.i,... the monic orthogonal poly- 



nomials with respect to \7.^ ) and {r.i}i=o,i.... the corresponding normalizations we can write 

p(a;o,a:2, . . . ,x„_2) = — TT (2rfe)qdet(r(x2j,a;2fc)) (7.7) 

O V / 0<j.k<n/2 

k—Q 



where 



T{x,y) := V 77—(xk{y)Xk{x)) = ( ^J'^'^] 

2rfe V J y D{x,y) S{y,x) 

( f{x)R2k{x) f{x)R2k+l{x) \ 

V !rfm2k{t)dt c fm2k+i{t) ) 

N/2-1 



Xkix) 



S{x,y) 



^ ^(i?2fc(2/) /" f{t)R2k+l{t)dt-R2k+l{y) I f{t)R2kit)dt 
I — n k J X J X 



k=0 
fV 

I{x,y) = I S{x,y')dy' 

J X 

D(x,y) = ^S^^^V) (7-8) 
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Proof. Because the polynomials {Rk{x)}k=o,i,... are monic we can add multiples of columns in to obtain 



p{xo,X2, . . . , Xn-2) = det 



f{x2j)Rk-l{x2j) 

rj{t)Rkit)dt 



n/2-1 

^ n (2^^) 



0<j<n/2 

0<fe<7l 



k=0 



X det 



/(a:2j)(2rfc)-i/2i?2^(a,2,) f{x2j){2rk)-^'^R2k+i{x2j) 



(7.9) 



0<j,fc<Ti/2 



Application of (7^) and the formula qdet A = qdet gives the formula (7/7) with S{x,y) as specified and 
formulas for I{x,y) and D{x,y) which are easily seen to be expressible in terms of S{x,y) as stated. □ 



A special feature of T(x,y), which follows from its definition in (7.^) in terms of Xk{y)Xk i^) ^^'^ skew 
orthogonality of {Rk{x)}k=o,i,... with respect to (7^), is the integration formulas 

T{x,x)dx = N/2 
T{x,y)T{y,z)dy 



T{x, z) 



(7.10) 



As a consequence of (7.10) and the quaternion formula (7.7) the integrations required to compute (B.l) can be 
carried out. Thus with ( 7.10 ) holding it is generally true that [|l6) 

/ c?a;2m qdet (T(x2i, a;2j) ) = (n/2 - (to - 1) ) qdet ( T(x2i, a;2j) ) (7.11) 

J — OQ ^ ^ 0<i,j<m \ / \ / 0<2j<m — 1 

Consequently we see from (^ that 

Pk{xo, ■ ■ ■ ,X2k-2) = qdet (T{x2i,X2j)) ■ (7.12) 

V / 0<ij<k 

If instead of considering even(OE„(/)) we consider odd(OE„(/)), the above working is essentially unchanged. 



Thus (7.7) and (7.8) hold with the replacements 

and 



{xq, X2,.-., Xn-2} {xi, Xs, . . . , Xn-l}, 

and with this modification of T{x,y) the formula ( 7.12| ) for pk holds with the replacements 

{Xo, X2,-.-, X2k-2} {Xl,---, X2k-l}- 



(7.13) 



(7.14) 



We remark that the structure of (7.12) with T{x, y) given by (7.8) is very similar to the general expression 
for Pk as computed for the ensemble SE„(((7//)^). First it is necessary to introduce monic skew orthogo- 
nal polynomials {Qk{x)}k=o,i,... and corresponding normalizations {qk}k=o,i,... with respect to the skew inner 
product 

dx {g{x) / f {x)Y {u{x)v' {x) — u{x)v'{x) 



{u\v)i := 
We then have Q (see also Q) 

Pk{xo, . . . ,Xfc_i) = qdet {Ti{xi,Xj)^ 



a<i,j<k 



(7.15) 
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where 



S4,{x,y) 
Di{x,y) 



I Si{x,y) h{x,y) 

\ Di{x,y) S^iy.x) 

N-l J., s 

Mr, 



E -^{Q2k{y)j^{f{x)Q2k+l{x))-Q2k+l{x)^[f{x)Q^k{x) 



dy 



2q, 

Si{x,y')dy' 



Si{x,y) 



(7.16) 



n odd 



The PDF for the distribution even(OE„(/)) with n odd is given by 



(n-3)/2 

^ n /(^20det 



1=0 



0<j<n 



As in (7.9) we can introduce the monic polynomials {-Rj(a;)}j=o,i,... to rewrite this as 



C 



det 



fix2i)Rj{x2'i.) 

Subtracting appropriate multiples of the last column from the columns 0,1, 
element of the column in the final row then gives 

f(x2i)Rj{x2i) \ 



f{t)R^^,{t)dt) 



det 



where 



Rj{x) := Rj{x) - 



XT dt I 0<,<(„-l)/2 

0<j<n 

C^f{t)R,{t)dt 



2 so as to eliminate the 



(7.17) 



/_°°^/(i)i?„_i(t)dt 



Rn-l{x). 



(7.18) 



The determinant in ( 7.17 ) is formally the same as that in (7^). Thus in the case n odd p{xq, X2i ■ ■ ■ , Xn-s) can 

n^n-l, Ri^Ri (7.19) 



be written as in (7/7) but with 



and C 1-^ C for some normalization C". 

Now we can check from the definition ( [7.1§| ) that for j = 1, . . . , rt — 1 the polynomials Rj-i satisfy the skew 
orthogonality property (7.3). This means that the integration formula (7.11) again applies in this modified 
setting and consequently the fc-point distribution is given by 



Pk 



{Xo, ■ ■ ■,X2k-2) = qdet (^f°'^'^{x2i,X2j)^ 



0<i,j<k 



(7.20) 
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where f°'^'^ is defined as in (7. J) but with the replacements (7.19). In the case of odd(OE„(/)) the replacements 
( 7.13 ) and ( 7.19 ) must be made in ( 7.12 ) and ( |7.8| ), and the replacement ( 7.14 ) made in ( 7.20| ). 



7.1 Summation formulas 



It has already been remarked that pk for SE„((g//))^) has the quaternion determinant form ( 7.15 ) and ( 7.16 ). 
Furthermore it is known Q that with / one of the classical weights in (^?l|), the quantity 5*4 in (7.16) can be 
summed to give an expression independent of the skew orthogonal polynomials associated with g, and dependent 
only on the monic orthogonal polynomials {pi(a:^)}i=o,i.... associated with the weight function g{x). Explicitly 



.9(a:)\i/2/(y) 



S2(x,y) 



ni—>-2n 



l2n-lf{y)P2n{y) / f {t)p2n-l{t) dt, 



(7.21) 



^9{y)^ fix) 

where 5*2 is specified by ( |6.8| ) and 72^-1 by Lemma 5.2, Here we will use results from ^ to obtain an analogous 



summation for the quantity S{x,y) in (7.S) 
Suppose n is even and write 



^ /"OO 

2 



^3 i^)-^^ I /(Osgn(x - t)Rj {t) dt. 
Then straightforward manipulation of the definition of S{x,y) allows it to be rewritten 

S{x,y) = '^{Si{x,y) - 5i(c»,y)) 



(7.22) 



where 



n/2-1 



Si{x,y)= V ^-^U2k[x)R2k+i{v)~^2k+i{x)R2k{y) 



fc=0 



The quantity Si{x,y) occurs in the quaternion determinant formula for fc-point distribution of 0E„(/). With 
/ one of the classical forms (pj) it can be summed to give |l| 



e / ^ f9{x)\^/^f{y) „ , . 
From this it follows 



+ 7n-2/(y)p«-i(y)^ 

ni— i-n— 1 Z 



sgn(a; - t)/(t)p„_2(t) dt. 



(7.23) 



S'i(oo,y) = 7„_2/(y)p«-i(y)2 / sgn{x -t)f{t)pn-2{t)dt 



and so by ( |7.22| ) we can evaluate S{x,y). 



Proposition 7.2. For {f,g) a classical pair \5.\ ), {pj(a;)}j=o.i,... monic orthogonal polynomials with respect 
to the weight function g{x), and n even the quantity S{x,y) in ^7.^ ) has the evaluation 



2S{X, y) = (^) "'lj^S2{x, y) - ln-2f{y)Pn-l{y) I f{t)pn-2{t) dt 



(7.24) 



(c.f. ^). 
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This summation fully determines even(OE„(/)) with n even. For odd(OE„(/)), n even, the prescription 



( 7.13 ) says the replacement i-^ should be made in ( 7.24 ). 

It remains to consider the case n odd. Consider first even(OE„(/)). In fact the formulas in ||l| giving the 
analogous formula to ( [7!23|) for n odd allows us to deduce that the summation ( [7.24[) remains valid for n odd. 
For n odd comparison of ( [7^ and shows 



S{x,y) = S^ix^y) 

which because of the formulas ( |7.8| ) (with nt-^ n — 1), ( |7l2| ), and ( |7l5| ), (fJ^ implies 

ovcn(OE2,i+i(/))/ X SE„((g//)^)/ x 



(7.25) 



(7.26) 



This is equivalent to the second statement of (4.87), which we already know from Theorem 5.1 is valid for the 
pairs (/,(?) in ( |5.l| ). In the case of odd(OE„(/)) with n odd, again the prescription ( [7.13 ) says we simply make 



the replacement I-oo '-"^ (''•2^) 
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